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Cosmology has provided a successful paradigm for 
structure formation in our universe through an inflation- 
ary phase [l| in early stages. Conceptually, however, the 
scenario is incomplete due to the presence of past singu- 
larities [2]. At such a singularity, the classical theory of 
general relativity breaks down and has to be replaced by 
an extended framework which remains well-defined even 
at very high curvatures. Since this requires modifications 
to general relativity at early stages of cosmic evolution, 
there can then also be corrections to the usual scenario 
of structure formation which might eventually be observ- 
able. While dimensional arguments and low energy effec- 
tive theory indicate that effects are very small, given by 
the tiny ratio of the Planck length lp = \J Gh to the Hub- 
ble length H~ x , a detailed analysis is required and may 
reveal more sizeable effects. This is what we provide in 
this letter in the framework of loop quantum gravity Q , 
a non-perturbative background independent approach to 
quantize gravity. 

Loop quantum gravity is one of the approaches where 
singularity resolution has been investigated using loop 
quantum cosmology [i| which results in the resolution 
of singularities in various situations including inhomo- 
geneous ones Semiclassical bounce pictures 
in special models have been described in 0, 0, Hi], [12] • 
A key role is played by the underlying quantum nature 
of spatial geometry [l3| ■ With such a discrete structure 
underlying classical space-time, effects not captured by 
low energy effective theory become possible. In partic- 
ular, there are large dimensionless parameters, such as 
the number of spatial lattice sites in a discrete state, 
which can always spoil dimensional arguments. In such 
a context, orders of magnitude of quantum corrections 
can only be estimated with a detailed analysis of the ef- 
fective equations arising from quantum gravity. Suitable 
techniques going beyond low energy effective theory are 
now available and are applied here. 



On larger scales farther away from the classical singu- 
larity one can use effective equations for the behavior of 
inhomogeneous perturbations, i.e. equations which are of 
classical type but amended by quantum correction terms 
as known from effective actions. Such equations can be 
used to test the semiclassical viability of crucial ingredi- 
ents of quantum gravity. At the same time, they allow a 
detailed study of the formation and evolution of structure 
including quantum corrections. 

In earlier papers , , candidate effective equations 
have been used for inhomogeneous scalar fields on an 
isotropic metric background. This allowed preliminary 
indications, but no systematic derivation or reliable pre- 
dictions. It therefore remained unclear which effects are 
to be expected. There are, for instance, cancellations 
in metric components on an isotropic background which 
hide terms that may or may not be modified by quan- 
tum gravity effects. Also, only corrections in the non- 
gravitational part of the field equations were considered, 
while the more complicated gravitational part plays a 
crucial role, too. Ignoring such corrections does not only 
affect details but may even change aspects such as the 
scale invariance of the expected perturbation spectrum 
which is already highly constrained observationally. 

For a reliable evaluation of this framework it is there- 
fore essential to derive full perturbation equations which 
take into account corrections of gravitational dynamics. 
This is reported here for scalar metric perturbations, 
which in Newtonian gauge are of diagonal form 

6ds 2 = ~2a{r 1 ) 2 <i>(r ] ,x)(dr] 2 +5 ab dx a dx b ) (1) 

on a flat background and in conformal time r\. Other 
gauges and modes can be included similarly, although 
we do not do this here for the sake of simplicity. Al- 
though quantum gravity is not formulated for classical 
metrics, the form of the metric (JIJ plays the role of select- 
ing the corresponding quantum regime where an effective 



2 



description is derived. This happens by picking semi- 
classical states of quantum gravity which are peaked on 
the given class of metrics, i.e. expectation values of met- 
ric operators are of the prescribed form and fluctuations 
around these values are small. Expectation values of the 
Hamiltonian operator in those states give the Hamilto- 
nian of the effective theory and thus effective equations 
XL |- Although the full quantization is background 
independent and non-perturbative, which is crucial for 
some properties of the quantum theory such as its spatial 
discreteness, a cosmological background is introduced in 
evaluating the theory through states and effective equa- 
tions. This puts the scenario in the usual context of 
cosmological perturbation theory, albeit including quan- 
tum corrections. Being derived from general semiclassical 
states which are not Lorentz invariant unless one restricts 
oneself to a vacuum state, effective equations may not be 
manifestly covariant even if the underlying quantum the- 
ory is covariant [l7j |. 

Compared to the standard derivation of cosmological 
perturbation equations [T3|, loop quantum gravity is in 
a different situation because it is based on a canonical 
quantization. Lagrangians are thus not quantized di- 
rectly, but Hamiltonians are used which provides an al- 
ternative but fully equivalent classical formulation. At 
this step, no new physics enters but it makes the formu- 
lation of quantum gravity possible. To apply this to the 
question of interest here requires a derivation of perturba- 
tion equations in a canonical framework, which allows one 
to include effective quantum modifications which arise for 
the effective Hamiltonians. Correction terms in the evo- 
lution equations then follow uniquely, but in an indirect 
manner. Thus, basic quantum modifications can have 
several complicated effects in the perturbation equations 
which allow one to test the underlying theory in a non- 
trivial manner. 

For loop quantum gravity in particular, quantum 
Hamiltonians are lattice operators taking into account 
the spatial discreteness of quantum geometry 
States are supported on lattices in space, which we as- 
sume here to be regular to simplify calculations; other- 
wise some coefficients can change but not by orders of 
magnitude. Lattice links are labeled by quantum num- 
bers p v j correspoonding to elementary areas (centered 
at a vertex v with transversal direction /) building up 
space. As seen from our final equations, the lattice does 
not introduce a preferred direction because the quantum 
Hamiltonian acting on the corresponding state is direc- 
tion independent. A Hamiltonian is constructed from ba- 
sic operators which are the elementary areas with eigen- 
values p v j and shift operators in those labels (related to 
curvature) [2l|. Variables are thus discrete, associated 
with lattices, and a classical geometry arises only in a 
continuum limit. The size of the p v j is given by the 
state as a multiple of £p, which is a general parameter 
constructed from G and H without input from quantum 



gravity. Minimum values where non-perturbative quan- 
tum effects are significant aiep v j ~ £p, but actual values 
of p v j in a semiclassical state can well be larger giv- 
ing rise to smaller quantum effects of the order £p/p v ,i- 
These are the terms whose cosmological implications we 
will study here. 

From lattice operators one first obtains an effective 
Hamiltonian, by taking an expectation value in a semi- 
classical state, as a function of the lattice areas rather 
than of the spatial metric. Coming from a lattice, such 
a function does not coincide with the classical one but 
contains discretization and other effects (whose magni- 
tudes are close to the extrinsic curvature scale) in addi- 
tion to £p/pv,i terms. Since p v j refers to the lattice size 
rather than the Hubble area, corrections can be much 
larger than £pH 2 as it was expected without discrete- 
ness. For p v j ss £p even non-perturbative quantum ef- 
fects have to be included, but for a semiclassical geome- 
try this cannot arise. On the other hand, those correc- 
tions become arbitrarily small for p v j — > oo, but there 
is an upper limit for p v j because large p v j imply large 
lattice sites. On length scales of ^/p^j, discreteness is 
noticeable which must thus be much smaller than scales 
probed by particle physics. Thus, a conservative upper 
bound is £p/y/p v ,i fCP 15 . But during inflation energy 
densities are much higher, up to Gp « f0~ 6 in Planck 
units, which requires £p/p v ,i 3> fCP 6 . This also means 
that extrinsic curvature terms given, through the Fried- 
mann equation, by */ Gp v jp and thus higher curvature 
corrections are small. Although precise estimates of cor- 
rection terms require detailed constructions of semiclas- 
sical states, the interplay of different corrections already 
suffices for a rough estimate of orders of magnitude. 

For regular lattices, all types of corrections can be de- 
termined explicitly for a given Hamiltonian [22^ . The 
Hamiltonian itself, however, is not fixed uniquely so far 
but subject to ambiguities such as the ordering of op- 
erators which do not commute in a quantum theory. A 
choice of Hamiltonian, including several ambiguity pa- 
rameters, corresponds to a fixed theory which can be 
tested phenomenologically. One can test precise aspects 
to constrain such parameters, or consider the whole class 
of possibilities allowed in the framework of loop quantum 
gravity and check whether this general behavior is viable 
at all. At the current stage, the second possibility is more 
reasonable to pursue and already very instructive due to 
tight constraints on the general properties of operators. 
This gives rise to the indicated quantum corrections when 
the resulting effective equations are expanded: First, spa- 
tial discreteness implies the replacement of differential 
by difference operators which, when expanded semiclas- 
sically on a background, result in higher derivative and 
higher curvature corrections. Secondly, inverse powers of 
metric components occur in Hamiltonians which would 
classically diverge near a singularity but are modified at 
small scales by quantum effects [23[. The latter correc- 
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tions of perturbative form t% /p v .i are most relevant for 
sub-Planckian curvature which is our present focus. 

Correction functions in coefficients of the effective 
Hamiltonian thus depend on the lattice areas p v j rather 
than a continuous field such as the spatial metric. More- 
over, they depend on time only implicitly through the 
time dependence of p v j. Qualitatively, such a correction 
function a(p v j) behaves in a way which approaches clas- 
sical behavior a = 1 for large p„,z but leads to suppres- 
sions of otherwise diverging inverse powers for small p v j 
241] . Most important for us is that any correction func- 



tion increases for very small p v j, reaches a peak of height 
larger than one and then approaches the classical expec- 
tation a = 1 from above in a perturbative expansion in 
^p/Pv.i- For perturbations around an isotropic geometry, 
one can express these corrections as functions of H since 
p~\ w Af 2 ^ 3 H 2 for TV lattice sites of volume p 5 Jf in a 
Hubble volume H~ 3 . The large factor A/" 2 / 3 thus magni- 
fies all corrections £pH 2 expected in low energy effective 
theory. From the perturbative metric (TTJ) it follows, on 
the other hand, that N 2 l 3 p v j{rj) cx a(?7) 2 (l - 2$(rj > v)) 
which allows one to write all effective equations in terms 
of the scalar perturbation <£>. 

The relevant gravitational dynamics is determined by 
the Hamiltonian 
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quantum gravity. Here, Ef is related to the spatial met- 
ric q a b by EfE\ — q ab det q and K % a is the canonical mo- 
mentum of Eg (related to extrinsic curvature). The con- 
nection r* depends on spatial derivatives of Eg and its 
inverse. The lapse function N is a free function but will 
be specified when choosing a gauge. 

For the effective Hamiltonian we keep only corrections 
for inverse powers of metric components, disregarding 
higher curvature corrections as it is adequate for sub- 
Planckian densities. There are two contributions by in- 
verse powers of the fields, the one explicit in |2]) and the 
other in connection components. We thus have two cor- 
rection functions, a multiplying the whole integrand and 
[3 multiplying connection components. They appear in 
different terms and will play quite different roles. There 
are thus classes of correction functions whose structure is 
determined theoretically and whose parameters, describ- 
ing their precise shape, can be restricted by observations 
or other means such as internal consistency. Despite of 
the non-uniqueness in parameters, crucial modifications 
are thus characteristic in the general form. 

As in any Hamiltonian system, the Hamiltonian gener- 
ates equations of motion. For the corrected Hamiltonian 
with N = a(l — $), corrected perturbation equations of 
scalar modes in conformal time rj take the form [25| 



a 2 /3V 2 $-3H$-3 1 - 



a p 



nH = --apST^, (3) 



expressed in basic fields (Eg , K 3 b ) which occur in loop 
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9 Q ($ + - 2a'p/a)) = ~f P<5T a ° . (5) 

A prime denotes derivatives with respect to p = a 2 which 
is the sum of all p v ,i, H = a/a, k = 8wG, and 5T£ are 
perturbations of the stress-energy tensor. For classical 
values of the correction functions, a = 3 = 1, we obtain 
the classical perturbation equations [l8(, which demon- 
strates the correct classical limit on very large scales of 
the effective theory. On intermediate scales, however, 
there are corrections which may lead to detailed and reli- 
able viability test of cosmological scenarios in loop quan- 
tum gravity and in proposals for potentially observable 
effects. 

Our effective equations include gravitational correc- 
tions from quantum gravity directly related to its ba- 
sic discrete structures. They suggest several applications 
on different scales, perhaps allowing tests of different 



regimes of quantum gravity. First, we isolate gravita- 
tional effects by assuming an effective perfect fluid back- 
ground such that T aa = P + SP — wTqq — w(p + 6p) with 
a constant w and T$ a — (p + P)u a with the energy den- 
sity p, pressure P and velocity u a . For curl- free velocity 
u a = d a u, we combine (|3]) and ([5]) to a Poisson equation 

V 2 $ - 5fi(p) 2 <S> = w^(Sp + 3a -1 W(p + P)u) (6) 
Zap 

with 8fi(p) 2 = 3H. 2 a'p/a 3 f3. For classical values, dfi = 
and we obtain the general relativistic Poisson equation 
corrected only by a pressure term. The correct classical 
limit is thus obtained as a — > 1 for £p — > 0. But since £p 
is non-zero, quantum effects always remain: to leading 
order we derive Newton's potential which recieves quan- 
tum corrections. Our derivation of the Newton potential 
is alternative to that proposed in [27| and conceptually 
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quite different. The precise value of corrections depends 
on the Hubble parameter, or the cosmological constant. 
This refers only to perturbations around a flat isotropic 
cosmology as this is the setting in which we derived our 
equations. The result on the Newton potential thus does 
not directly apply to the solar system for which perturba- 
tion equations around the Schwarzschild solution are re- 
quired. They can be derived by the same methods which 
are, however, technically more involved for a curved and 
inhomogeneous background. From the general procedure 
we expect that the Hubble parameter occuring in correc- 
tion terms will effectively be replaced by the solar mass. 

Additional applications arise for structure formation. 
For the main effect we combine Eqs. (j3|) and (J4j) to elim- 
inate stress-energy, still for an effective perfect fluid: 

$ + 3(1 + w + ei)W$ -{w + e 2 )V 2 $ + e 3 H 2 <P = (7) 

with quantum corrections e^. The equation is sensitive 
to the gravitational part of perturbation equations whose 
correction terms, derived here for the first time, require 
quantum gravity. This equation has a characteristic im- 
plication: classically, = and the last term cancels 
exactly, but with quantum corrections e 3 = —2a"p 2 /a is 
negative. These corrections are usually small compared 
to the term V 2 <£>, but can become important for modes 
of small comoving wave number k such that ak <C H. 
Such modes are outside the Hubble radius for which 
classically <£> would be preserved. With quantum cor- 
rections, however, such curvature perturbations are no 
longer conserved. This is seen by solving the equation 
$ + (l + ^)<j>/r/ + e 3 $/?7 2 = 0, with v = (5 + 3«j)/(l + 3to), 
for the behavior of large scale modes in conformal time 77, 
approximating w and e 3 by a constant for the qualitative 
effect: $(77) = 77 A with A = -f ± \^/v 2 - 4e 3 . Clas- 
sically, there is one decaying mode and a constant one, 
corresponding to conserved With non-zero e 3 , how- 
ever, the constant mode disappears, affecting the power 
spectrum. Since not only the magnitude and possibly 
the sign of e 3 but also cosmic evolution time depends 
on the mode, a specific running of the spectral index 
can be expected. We comment here only on the mag- 
nitude of corrections for large scale modes which were 
created early in inflation: as before, an estimate for a—1 
gives 1 ^> |e 3 | ^> 10 -6 . During inflation, conformal time 
77 oc e~ Hr for modes currently visible on the largest scales 
changes by a factor e -60 , such that the constant classical 
solution is corrected by a factor e _60e3 «1- 10 2 e 3 . On 
the lower end of e 3 this would not be observable soon, 
but with lattice areas expected to be much closer to the 
Planck scale, e 3 should be closer to one and the magnifi- 
cation due to the number of e-foldings further enhances 
quantum corrections to become potentially observable. 

We have provided first effective cosmological perturba- 
tion equations which include correction terms from quan- 
tum gravity. The scheme of the derivation is system- 
atic and general enough to include also other modes and 



backgrounds. Unlike low energy effective theory which 
is usually used to introduce quantum corrections in clas- 
sical equations, effective theory taking into account the 
spatial discreteness expected from quantum gravity has 
revealed new effects whose magnitude can be much larger 
than expected on dimensional grounds. In addition, ef- 
fects can be enlarged due to long cosmic evolution. How 
large precisely corrections will be requires more detailed 
solutions of lattice states, possibly involving new numer- 
ical schemes, which is accessible in the framework now 
provided by loop quantum gravity. 
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